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Thermal fluctuations cause the local normal vectors of fluid interfaces to deviate from the vertical 
, direction defined by the flat mean interface position. This leads to a nonzero mean value of the 

. corresponding polar tilt angle which renders a characterization of the thermal state of an interface. 

Based on the concept of an effective interface Hamiltonian we determine the variances of the local 
interface position and of its lateral derivatives. This leads to the probability distribution functions 
for the metric of the interface and for the tilt angle which allows us to calculate its mean value 
and its mean square deviation. We compare the temperature dependences of these quantities as 
predicted by the simple capillary wave model, by an improved phenomenological model, and by the 
microscopic effective interface Hamiltonian derived from density functional theory. The mean tilt 
angle discriminates clearly between these theoretical approaches and emphasizes the importance of 
the variation of the surface tension at small wave lengths. Also the tilt angle two-point correlation 
f j — ^ , function is determined which renders an additional structural characterization of interfacial fluctu- 

Q ■ ations. Various experimental accesses to measure the local orientational fiuctuations are discussed. 



PACS numbers: 68.03.-g, 68.37.-d, 68.05.Cf, 68.15.-|-e, 68.47.Pe 



INTRODUCTION 



Q If two phases of condensed matter, which coexist thermodynamically at first-order phase transitions, are brought 

^ into spatial contact under suitable boundary conditions, an interface is formed through which the system interpolates 
smoothly between the structural properties of the adjacent bulk phases. The formation of such an interface is 
^ ■ accompanied by a cost in free energy given by the product of the mean interfacial area and the surface tension. 
^ ] Since the surface tension is always positive the system acts as to minimize the interfacial area under given boundary 
■T-H ■ conditions which gives rise to a plethora of capillary phenomena (for reviews see, e.g., Refs. HSiSi). If the 
^\ ' first-order phase transitions end at a critical point where the differences between the coexisting phases disappear, the 
interface broadens and ceases to exist; accordingly the surface tension vanishes there with a universal power law as 
lO function of temperature. There are numerous experimental and theoretical techniques to measure and to calculate 
J— ^ , surface tensions the knowledge of which allows one to tune capillary phenomena in a specific way 0,1111 Si. 

With the advent of sophisticated spectroscopic techniques such as X-ray and neutron scattering under grazing 
incidence it has become possible to resolve the structural properties of interfaces at an atomic scale, paving the way 
"j^ for a microscopic understanding of interfacial properties in terms of mean density profiles, concentration profiles, and 
- in the case of nonspherical particles - orientational profiles (see, e.g., Ref. 'Gf and references therein). 
I These mean profiles, i.e., one-point correlation functions, can be inferred from deviations of reflectivity coefficients 

'■^ ' from the corresponding Fresnel formulae holding for steplike variations of the structural properties. The analysis 
Ch . of the diffuse scattering around the specular beam provides an even deeper insight into the structural properties of 
O I interfaces by probing the structure factor and thus the two-point correlation function in the lateral directions and 
close to the interface 

This interfacial structure factor is of particular interest because in addition to the bulk-like density fiuctuations, 
which occur on length scales up to the bulk correlation length ^ and reach the surface of each of the two semi-infinite 
phases, there are capillary wavelike fiuctuations occuring at length scales larger than ^ but tied to the interfacial region 
only. Without gravity or lateral confinement these capillary wavelike fiuctuations lead to a unlimited thickening of 
the interface and to a unlimited increase in the spatial dependence of the two-point correlation function on the lateral 
distance, known as the roughening of the interface. If one of the two adjacent phases is a solid, this broadening 
occurs only above the so-called roughening transition whereas for interfaces between fluid phases it occurs at all 
temperatures. Thus, in the following we focus on fluid interfaces, in particular the liquid-vapor interface of a simple 
fluid. 

The physical picture of the capillary wavelike fluctuations is that of an intrinsic proflle - generated by bulklike 
fiuctuations as if the interface would be fiat - exhibiting local, laterally varying height positions. This picture reconciles 
van der Waals' original idea of an interface |3and the purely capillary wavelike picture of a fluctuating steplike interface 
put forward by Buff, Lovett, and Stilhnger Q. 

Based on microscopic density functional theory one can determine the effective interface Hamiltonian, and thus the 
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statistical weight for such a given configuration, which can be expressed in terms of a wavelength dependent surface 
tension; it equals the macroscopic surface tension at large wavelengths, attains a deep minimum at medium lengths, 
and increases for short wavelengths . The ensuing theoretical predictions for the lateral interfacial structure factor 
have been confirmed quantitatively for water [3, OMCTS, CCI4, squalane, glycol ^.nd Ga by using X-ray 
scattering under grazing incidence. 

Given this successful agreement between theory and experiment one can turn to the question whether this wavelength 
dependent surface tension represents primarily a transparent parametrization of the interfacial structure factor or 
whether the corresponding effective interface Hamiltonian has the quality of transferability, i.e., predicts also other 
key features of fluid interfaces. 

In this context we consider a snapshot of an interfacial configuration. The points in space, at which the local 
density equals, say, the mean value of the two bulk densities, form a two-dimensional isodensity surface which can be 
identified as the local interface position. This surface is characterized by giving its local and laterally varying height 
with respect to the mean flat interface position. 

Alternatively, at each point of the surface one can attach its normal vector which forms the polar angle 6 with the 
vertical direction given, e.g., by the direction of gravity and which is the normal for the flat mean interface. The 
thermal fluctuations of the interface cause corresponding fluctuations of the surface normals and thus of the polar tilt 
angle 9. Since 9 is always positive its mean value < 9 > is nonzero; < 9 > — would imply the complete absence 
of capillary wavelikc fluctuations. Therefore, < 9 > and the mean squared deviation < 9^ > — < 9 >^ characterize 
the thermal state of the interface with a particular emphasis on the capillary wavelike fluctuations. Based on the 
aforementioned effective interface Hamiltonian one can derive the probability distribution function for 9 and thus 
calculate its mean value and its moments. Their increase with temperature T signals the unfreezing of capillary 
waves. In this sense < 9 > (T) represents an additional characteristic quantity for defining the thermal state of 
interfaces, on equal footing with the surface tension. Moreover, the function < 9 > (T) probes directly the concept 
of the wavelength dependent surface energy. 

Unrelated to this latter aspect, in a series of papers Simpson and Rowlen studied the surface roughness contributions 
to measurements, based on linear dichroism, of the molecular orientations of optically active molecules adsorbed on 
rough solid surfaces [T3LIT3.IT5l| . In order to infer from such data the local molecular orientation they must be corrected 
by the effects due to the frozen oricntational fluctuations of the rough supporting substrate. Thus in this context the 
tilt angle distribution of the interface plays an important role. Other experiments consider oricntational fluctuations 
in monolayers or thin multilayer films of non-spherical particles adsorbed on fluid interfaces jl6, 17, 18, 19, 2(1] . 
These oricntational fluctuations are determined by fixed relations between the orientations of the adsorbed molecules 
and the fluctuating local orientations of the supporting fluid interface. Here we focus on the latter contribution. 
There is a wide range of experimental techniques which provide access to the molecular orientations at interfaces [T^ 
such as fluorescence, second-harmonic generation (SHG), sum- frequency generation (SFG), linear dichroism (LD), 
and angle-resolved photoacoustic spectroscopy (ARPAS). For instance, in Ref. the spectroscopically determined 
orientation of complex molecules such as proteins attached to rough solid substrates is studied. If adsorbed with low 
concentrations at fluid interfaces such tracer molecules probe the oricntational fluctuations of the interface normals. 
Promising further perspectives are provided by studying fluid interfaces of colloidal suspensions by laser scanning 
confocal microscopy which allows one to observe, down to the scale of the colloidal particles, thermally induced 
capillary waves in real space and to analyze the distribution of the local interface orientations '22] . 




In Sec. njwe introduce the geometry of the local orientations of fluid interfaces followed in Sec. Illll bv a discussion 
of the various theoretical approaches describing capillary wavelike fluctuations as they are used for calculating the 
Gaussian interfacial fluctuations in Sec. IIVI The dependence of the mean tilt angle and its mean squared deviation 
on physical parameters is analyzed in Sec. [3 The relation of these results with experiments is discussed in Sec. IVII 
We summarize our results in Sec. IVIII 



The capillary wave approach to study liquid interfaces, which was put forward by Buff, Lovett, and Stillinger (BLS) 
, describes the actual smooth density profile as the thermal average of a locally fiuctuating steplike interface between 
the two coexisting phases. The local position of the liquid-vapor interface is described by a single-valued function 
/(R) depending on the lateral coordinates R — {x, y) with R — |R| in the x-y-plane parallel to the mean interface at 
z = (see Fig. 1). This so-called Monge parametrization neglects overhangs of the interface as well as bubbles of one 
phase inside the other, i.e., domains topologically separated from the interface; apart from the close neighborhood of 
the critical point such configurations carry a sufficiently small statistical weight to be negligible. The orientation of 
the mean interface described by the unit-vector is chosen to be parallel to the z-axis, i.e., perpendicular to the x-y- 
plane. This direction is determined either by corresponding boundary conditions or the direction of the gravitational 
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force. Although the mean orientation is fixed one does observe local variations of the interfacial orientation due to 
thermal fluctuation of /(R). These deviation may be described by a local tilt angle 



(R) = arccos(ez • n(R)) 



where 



n(R) 



-fy 



1 + + n \ 1 



(1) 



(2) 



denotes the normal vector of the interface at the lateral position R which is given by the partial derivatives fi := 
i = x,y, of the interface position /(R). In addition to the polar angle < 6'(R) < tt/2 one may define 



0(R) = arccos 



n(R) ■ e^: 
sin6i(R) 



(3) 



i.e., the azimuthal angle between the projection of the normal vector n onto the x-y-plane and the x-direction in the 
x-y-plane. The local tilt angle 6'(R) and 0(R) describe uniquely the interface orientation because 



^cos(/)(R) sin6'(R)' 
n(R) = ( sin(/)(R)sin6'(R) 
cos6'(R) 



can be reconstructed from the knowledge of the angles and 



(4) 



z = 



vapor (pg) 




G 



FIG. 1: Schematic picture of an interface configuration between the coexisting liquid and vapor bulk phases with number 
densities pi and Pg, respectively. The interface does not contain overhangs or bubbles. Thus the local position of the liquid- vapor 
interface can be described by a single-valued function /(R), where R = {Rx,Ry) denotes the lateral coordinates. The dashed 
curves indicate the intersections between the two-dimensional manifold /(R) and planes x = = const and y = Ry — const, 
respectively. Gravity G leads to a mean interface position at 2: = (dotted lines). The volume of the sample is V = AL where 
A is the lateral area and L is the vertical height. The local orientation of the interface is described by the normal vector n(R) 
characterized by the angles 0(R) and 0(R) (see Eq. 



4 



An important motivation to study the local tilt angle < 0(IV) < tt/2 is the non-vanishing thermal average 

<0> (T) > (5) 

of the tilt angle which is strictly positive at finite temperatures T. Thus the interface is always locally tilted, so that 
one may observe relevant thermal fluctuation features already in the mean value < 9 > (T), in contrast to the mean 
interface position < / > (R) = which always vanishes. A vanishing mean value of 6 would imply the complete 
absence of capillary wavelike fluctuations. Thus, the 'structure' of the interface reveals itself via the interface position 
/(R) only through its variance, i.e., the mean squared deviation 

a\T) = <(/(R)f > = J^^^^a^iq) (6) 

from its vanishing mean value, where o'^iq) denotes the fluctuation contributions stemming from interfacial modes 
with lateral wavevector q with q = |q| (see, c.f., Eq. 123). 

Recent progress in grazing incidence X-ray scattering experiments [Tol llll Il2l | enables an almost complete deter- 
mination of (7^(q) and therefore of the structure and the fluctuations of the interfacial position /(R). In particular, 
as predicted theoretically , it has been demonstrated by X-ray scattering experiments that for wavevectors above a 
few nm^^ there is a decrease of up to 75% of the surface energy due to the attractive dispersion forces. But for some 
fluids the fluctuation spectrum at large wavevectors q could not be resolved close to the physical limit qmax = v 
a few where a denotes the microscopic size of the molecules. As we shall show in the following the fluctuations 
of the local tilt angle 6'(R), i.e., its mean value < 6 > and its mean squared deviation < {69)'^ > with 69 — 9~ < 9 > 
are much more sensitive to the structure of the interface at large wavevectors, so that the analysis of them promises 
a significant improvement of the understanding of fluid interfaces. 



III. CAPILLARY WAVE THEORY 



The cost in free energy to generate an undulation of a flat interface, i.e., to generate a non-flat interface configuration 
depends on the amplitude /(R) of the undulation but also on physical parameters such as temperature T, the 
gravitational constant G, the surface tension 70, and the bending rigidity kq or ultimately the interaction potential 
between the fluid particles. The usual approach to describe this cost in free energy assumes that the amplitudes are 
small so that the Hamiltonian may be expanded into powers of /(R) and its derivatives V/(R) and A/(R) yielding 
(m denotes the mass of the particles and Ap = pi — pg the difference in number densities of the coexisting liquid and 
vapor phase) 

7^(^)[/(R)] SR LgAp(/(R))^ 

^ JR2 L (7) 

+70 (V/(R))' + «o (A/(R))' . 

This gradient expansion is known also as the Helfrich Hamiltonian for fluctuating membranes |U 13113 ■ 

It is transparent to study 7i^^' in Fourier space in which the bending modes decouple. To this end we introduce 
the Fourier transformed functions 

/(q) = / d?Re-^^-^f{^) , 

(8) 

/(R) = / A^^e^^-^m 
yielding the Fourier transformed Hamiltonian 

n^"\m] - (0 {r^^GAp + ^^^Hlh') (9) 

with the momentum dependent surface tension 

I'^^Hq) = 10 + lioq^ ■ (10) 

We note that the use of an upper cut-off q < Qmax = ^w/a for the wavevectors is indispensible because the concept of 
capillary wave fluctuations becomes meaningless on length scales smaller than the microscopic size a of the molecules. 
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Instead of starting from the ansatz in Eq. (O an alternative strategy (see Ref. 0) to derive 7i[/(R)] is based 
on a microscopic density functional theory for inhomogeneous simple fluids, which is a successful approach for the 
description of nonuniform fluids [2^ f27l | . The different kinds of occuring density fluctuations - bulk bubbles and 
interface undulations - arc separated by determining the intrinsic density profile Pf{r) via minimizing the grand 
canonical density functional r2[/9(r)] under the constraint of a locally prescribed interface position /(R), i.e., the 
location of the isodensity contour of the mean density is fixed as function of the lateral coordinates R. Due to 
the long-ranged character of the underlying van der Waals interparticle potential W{r) between the fluid particles 
nonanalytic contributions occur and therefore a gradient expansion of the Hamiltonian breaks down. Instead, in 
Fourier space one obtains in terms of the amplitudes /(q) describing the interface position the Hamiltonian 

n^'^Hm] - I S^,\E{q)\m? (11) 

7r2 (27r)^ 2 

where the superscript (G) indicates that this expression refers to a Gaussian approximation bilinear in the amplitudes 
/(q) of the bending modes. The momentum dependent energy (per (lenght)*) 

E{q) ^Eo + l[q)q^ (12) 

with Eq = mGAp can be written in terms of a momentum dependent surface tension 'y{q), which is not simply given 
by the Helfrich expression H10|l but by an explicit functional of the molecular interaction potential VF(r). Based on 
the simple functional |2^ 

n [pir]] = [ d\ ihipiv)) + pp{v) + p{v)V{v)) 
Jv 



+ i / d\ I d\'w{\v-v'\)p{v)p{v') , (13) 



V Jv 



where V is the volume of the sample, p(r) is the number density of the fluid particles at r = {x,y,z), r = |r|, fhip) 
is the reference bulk free energy of a system determined by the short-ranged, repulsive contribution Ws{r) to the 
interaction potential W(r), and w{r) is the attractive part of W{r), one finds within the reliable so-called product 
approximation 



^w"{0)q^ 



with the surface tension and bending rigidities 



70 = ^^"(0)^ 

4""^ = ijw"io)iApr^cH > 



24 



«o = ^(Ap) e Ch^-^ > 
For the bulk free energy density fh{p) one may use the Carnahan-Starling expression [27 



(14) 



(15) 



Mp) = ksTp |ln(pA3) - 1 + ^^3^} (16) 

where A is the thermal de Broglie wavelength and rj — ^prQ the packing fraction. (For the definition of tq see, c.f., 
Eq. (|18|l ). Within this product approximation the momentum dependence of 7(g) in Eq. 114|l is determined by the 
three-dimensional Fourier transform 

wiQ) = / d^r e-'*5 ""w(|r|) , (17) 

JR3 
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of the interparticle potential w{r) evaluated for the absolute value q of the lateral momentum q with w"{0) > and 
w{0) < 0. 

Considering particles which interact via van der Waals forces one may adopt for the attractive part of the interaction 
potential w{r) the form 



w{r) 



(r2 +r2)3 



(18) 



which exhibits the correct large distance behavior w{r — > oo) — > —Ar^^ with A = woTq for c? = 3, r = 3 28]. The 
length ro corresponds to the diameter of the particles and thus serves as a lower limit for the length scale of the 
density fluctuations and of the capillary waves considered below. Then, Eq. H14() reduces to the explicit expression 



70 [qro) 




0.74C|,(Cg)'U+(^) (1 - (1 + gro)e-''-«) 



3 V ^0 



(19) 



0.74C|f 




{qro)' - Ch 



{qro 



O (q^e- 



depending solely on the bulk correlation length ^, the diameter tq of the molecules, and the dimensionless parameter 
Ch ~ 0.25 describing the distortion of the intrinsic density profile due to the bending of the interface. Upon increasing 
Qi i{q)/1o decreases below 1, reaches a minimum and increases again (for more details see Ref. 0). This decrease of 
7(g) is in accordance with recent simulation data |2^|33 which, however, have not yet confirmed the predicted and 
experimentally observed reincrease of 7(9) at large q. 

For the density functional theory given in Eq. (|13l) and the interaction potential in Eq. IjlSf) one has w{q) = 
■«;(0)(1 + qro)e^'^^° with w{0) = — 7r-^rQU)o/4, so that one finally finds the temperature dependences 



70 



48 



wof{Ap)' 



~{H) 

Kg 



(20) 



and close to the critical point at T^: 



f.iHH)l Irl 



Ap 



1.23 



T 



(21) 



2Tr\ Tr 



Although it is straightforward to determine ^ and Ap numerically we found that these approximate but analytic 
expressions are quite accurate for temperatures between the triple point Ttr — {2/'i)Tc and the critical point (apart 
from the critical region in which deviations from the above mean- field behavior occur), so that they can be used in 
the following. For temperatures below Ttr - where no fluid interface exists - the model defined by Eq. H13|) ceases to 
be applicable because it does not capture the freezing transition, which would require a more sophisticated version of 
the density functional. 

In order to study the dependences on the numerous physical parameters we define two length scales: the capillary 
length Lc and the stiffness length L^, respectively. 
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which in addition to the microscopic cut-ofF lengths tq (range of the hard-core repulsion potential) and a = ^'^ (size 
of the molecules), and the bulk correlation length ^ are relevant for the cost in free energy for bending an interface 
with wave vector q (see Eqs. (Illll and (|12|l ). For simple liquids the values for ro, a, and Lk are of the same order of 
magnitude and may be identified in a first approximation. But since the distributions of the local orientations of the 
interface turn out to depend sensitively on these microscopic length scales we prefer to keep them distinguished until 
their influence has been elucidated. Moreover, colloidal suspensions (see, e.g., Ref. jl^) might offer the possibility to 
keep these length scales well separated and to vary them to large extent independently. 



IV. GAUSSIAN INTERFACIAL FLUCTUATIONS 



Before we study in detail the dependence of the interfacial structure on physical parameters such as the surface 
tension 70 or the interparticle potential w(y) in this section we discuss those properties of an interface which follow 
from the Gaussian approximation given by Eq. 

Since within the Gaussian approximation P[/] ~ exp(— /3H^*^^), see Eq. (tTTll . the bending modes are decoupled one 
finds for the probability of a surface wave /(g) of wavevector q the explicit expression (see Eq. \\\\ \ 

nm - ^e-mm.)f (23) 

so that mean values < A > for functionals A[/(R)] of the interface position can straightforwardly be calculated by 

<A> [pE{q)] ^ I'D! ^[/(R)] (24) 

where 

fvf...= n (^^^~(^)) r ' H^'^^) ^e-mi-^m-H^-fi^m (25) 

denotes the integration measure for the interfacial degrees of freedom which solely depends on f3E{q) with /3 — 1/ (ksT) 
as the inverse thermal energy. Because of the relation /(— q) = /*(q), the product in Eq. H25|l runs only over suitably 
discretized vectors q in the half space M^/2. Details on the Wiener measure in terms of Fourier coefhcients can be 
found in Ref. 0. 

In a first step we consider local second order moments cr|„ 2m (see, c.f., Eq. H29|l l characterizing the distribution 
functions of lateral derivatives of /(R) at one and the same lateral position R, whereas Subsec. IIV Bl focuses on 
correlations of /(R) at points R = Ri, R2 with non- vanishing distance vector a = R2 — Ri. 



A. Local variances 



The probability P[fo] to find the interface at the lateral coordinate R at a specified position /o = /(R = 0) is 
given by (see Eq. (jSJ) 



P[/o;.]^/^/.(/o-/^^|^/(.) 



f ^^^Pfo Yl r d(nef{<i)) r d(lmfici)] 

-^(g) -|^[(7^e/(q))" + (Tm/(q))"] -2»p(Ke/(q)) 



1 JQ 

=e 2^ 



V2^ 

with the fluctuation spectrum (see Eq. JH))), o^ifi) = kBT/E{q) and the variance 



2 / d^q ksT f ksT q 
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Integrating over all wavevectors requires special care since modes may be constrained by physical boundary conditions. 
As stated after Eq. (|1U|) Qmax — 27r/a denotes the maximum possible wavevector because fluctuations on length 
scales smaller than the microscopic size a of the molecules are not compatible with the concept of capillary waves. 
Accordingly, one may introduce a lower cut-off qmin — 27r/L due to the finite lateral size L of the interface. If L 
is larger than the capillary length Lc (see Eq. (|22|) ') one may safely set qmin = so that fluctuations with long 
wavelengths are suppressed only by gravity. 

Indeed, P[fo;a] is properly normalized so that c(foP[/o; cr] = 1. Accordingly, the probability ^'[/q"'™^; cr2n,2m] 
to find a prescribed derivative 



(28) 



R=0 



at R is given analogously to Eq. 1)26(1 with /o replaced by /q"'™'' but with the variance 



2 _ / d, q ^2n„2m _ , , 

"2n,2m — / (2Tr)'^ E{q) ~ '"^ "2n+2m+l V^^/ 



with 



^a'^ (30) 

and In^m = Jq'' ^ cos^" (^sin^™ 4> ^ ^ {Jn,m + Jm,n) for integers n, m G N with Jn^rn = 2 /J"^^ d0 cos^" sin^™ (j) = 
B{m + i,n+ i) (3.621.5 and 8.384.1 in Ref. [H) so that (n, m e N) 

/„„.^l"(" + ^)"(-t^^ ■ (31) 
TT r(n + m + 1) 

In particular, one finds /o,o — 1, Ii.o = Io,i = |j ^2,0 — ^0,2 = |, — |- Therefore the probability P[fo;o''] = 
-ff/d^'^^ f''']^[/o°'^''; foi' the slope (orientation) (q — V/(R = 0) of the interface at R = has the variance 
a''^ — (^20 — '^i/'^- Note that cr', 173, and fp are dimensionless. The probability P[Af;aH] of the Laplacian A/(R), 
which in lowest order in / is twice the mean curvature 

fxxi^ + fy) ~ '^fxfyfxy + fyyi^ + fx) 

2(1 + /.^ + /,^)^ 

has the variance tr^ = : 

2 _ f d^l ( 2 2\2 ^s?" _ f dq ^ksT 



B. Finite distance variances 



An alternative to the local slope fg as a measure for the orientation of an interface is provided by finite differences 
5f := /(R + a) — /(R) of the interface position. Analogous to Eq. H26|l . the joint probability of finding the heights 
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/o and /a of the interface at the lateral positions and a 7^ 0, respectively, reads 

p[/o, /a] - / w ^ (/o - / ^m) ^ (/a - / ^fi^y^-' 

00 00 



^.-|5Tf('^e/(q))'' + (2^'»/(q))1-2i(Ke/(q))(p+p' cos qa)+2j(lm/(q))p' sin qa 
00 00 



-00 —00 



using the formula 



2/ N f (Pq ksT 



gmax 



(34) 



J dxi exp |--a;iAija;j + Sia;^! = --^== exp {s^ A~^Sj } (35) 

for any symmetric and positive definite matrix (Aij) with the matrices of the variances (S'^S^^ = I) 

where ct^ is given by Eq. H27|l . The eigenvalues 

A±(a) =cr2_^o-2(a) > (37) 
of the matrices and A+A_S^^ are determined by Eq. (|27l) and by the correlation function 



(38) 



with a'^{a = 0) = and where Jo{x) denotes the Bessel function with Jo(a; — + 0) = 1 — + 0{x'^) and Jo{x 
00) = (sin X + cos a;) + ©(2;^^). Since E{q) increases at least as ~ for large q the latter property of Jo{x) allows 

one to set qmax = 00 for any finite a. For qmin = and qmax — 00 and E{q) given by Eq. (|10|) one finds for the 
correlation function defined in Eq. 



ln(Lc/a) , Lc » a , 



(39) 



2vr7o 

where ifo(a;) denotes the modified Bessel function. One recovers the distribution of the interface position (Eq. (|26|l 'l 

00 

P[fo;'j]= I dU P[/o,/a;S] (40) 

— 00 

by integrating over the positions at a. 
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Note that the conditional probabiHty 

P[/0,/a] _°^P 



I 2 



i-y- , /f/^"-2/o/^o-"(a)/^- 
*(a)-l l-cr'i(a)/cr* 



0} 



V2w2(l-a4(a)/a4) 



(41) 



is not a function of 6f ~ fo — fg^ alone, but depends on /o and /a separately. For the distribution of the difference 
5/(a) := /(R + a) — /(R), i.e., the probability that at given points R and R + a the height difference has a certain 
value <5/(a), one finds not the conditional probability, Eq. I|41|l . but the expression 



oo 

P[(5/(a)] = j df P[fo + f,U + f] 



df 



1 



exp ■ 



/0+/a 

2 



4A. 



2A+A- 



(42) 



1 (££(a))^ 



with 



a^(a) = 2A_ = 2(cr2 - (j2(a)) = 



1 



dq 



q{l - Joiqa)) 



(43) 



Due to the translational invariance in lateral directions this distribution for /(R + a) — /(R) is independent of R 
so that without loss of generality one can set R = 0. The finite difference variance (7g{a) is shown in Fig. El for the 
phenomenological capillary wave model with E{q)/EQ = 1 + L^q^(l + L\q^) (see Eqs. (|10ll and (|12|) ') and the effective 
length scales Lc = ^/jo/Eq and — a/ko/to given by Eq. . 



2.5 




FIG. 2: Finite difference variance a"|(a) (Eq. M^jl l normalized by the variance of the interface position (Eq. I|27|l 'l as 
function of a/L^ for various values of L^/L^ (Eq. i22l ) for tire phenomenological capillary wave model (Eqs. ^ and ([HI). 
The quadratic increase of (T|(a) in the limit a ^ (see Eq. I^U) is visible on the scale g/Lk only for small values of Lk/Lc. 
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In the limit |a| ^ oo one recovers the variance 2cr^ corresponding to two independently fluctuating interfaces. In 
the limit |a| ^ one finds 

a2(a) ^ 3? + 0{{q^axaf) (44) 

in accordance with the expectation that for small distances a the difference /(R + a) — /(R) is determined by the 
slope of /(R), so that its variance is given by the variance a' = (T2,o = o'3/V2 of the slope (see Eq. times the 

distance a. 

The comparision of Eqs. H43|) and H3U|) allows one to express the variance cr'^ of the local slope in terms of the finite 
distance variance cr^(a) for the height distribution at an effective lateral distance a{Li^qmaxTLcqmax)'- 

a'^al ^ aj{a{Li^q„iax,Lcqrnax)) , (45) 

where ct'^Oq is a local estimate of the variance ag{ao) at a finite distance oq. For small values of uq one may apply 
the expansion in Eq. (|44|l yielding a{L^qmaxT Lcqmax) = ciq, but for finite distances the effective lateral distance 
o-iL^qmax, Lcqrnax) IS determined implicitly by 

max nQmaxi ^ cQmax^Q_max y)) 

{Lcq-max)''^ + {LnqraaxYv'^ 



„ / 7 y V y'^KjHmaxtf I I ^ '-'v\'-^\^K^maxi -"-^c^max j^maxy j t 

= / — ^-2 , .,2,(T ^274 ■ (46) 



The solution of Eq. H4t)|) is plotted in Fig. 13 for several values of the capillary length Lcqmax ■ For sufficiently small 
distances oq one finds a{L^qmax, Lcqmax) — lo- In general, the numerical equivalence of the finite distance approach 
and the local slope estimate, i.e., of a and ao becomes better for increasing length scales Lc and L^. 

In order to obtain a good agreement of the local slope estimate cr'^a§ and the finite distance variance (T|(a) the 
distance ao should be smaller than and also smaller than 27: / q^ax- 

The finite distance approach may be used as an alternative to regularize the fluctuation spectrum, i.e., to suppress 
the infrared divergence caused by the large wavelength mode q ^ Q which leads to an infinite variance a for vanishing 
gravity, i.e., Eq — 0. Instead of introducing a minimum value qmin for the allowed wavevectors - which is necessary 
for local quantities such as the variance studied in Subsect. IIV Al (see Eq. (|27|) ') - one can consider correlations at 
finite distances a for which one can set qmin = 0. 

Also the ultraviolet divergence, which is caused by the diverging number of large momenta modes and yields 
an infinite variance a' — 0-3 /-\/2 of the slope (see Eq. - even for a non- vanishing bending rigidity kq, may be 

regularized by the finite distance approach. Instead of introducing a maximum value qmax for the allowed wavevectors 
- which is necessary for local quantities such as the variances studied in Subsect. IIV Al fsee Eq. i27\i ) - one can consider 
differences /(R + a) — /(R) at finite distances a for which one can use qmax = 00 if kq 7^ 0. The corresponding 
variance (Jg{a) of the difference Sf{a) remains finite for kq > whereas even in this case the variance ct'^ of the local 
slope diverges. 

But as one can see from Eq. (I44f) the finite difference variance <Jg{a) is given by the variance cr'^ of the local slope 
for small values of a, so that both approaches are even almost quantitatively equivalent for distances a < ^'^ or 
for large qmax- At least for variances significant differences cannot be observed. Therefore in the following we apply 
the simpler and conceptually more straightforward approach of local orientations with wavevectors larger than qmax 
being not permitted. 



C. Distribution of metric g and tilt angle Q 



The geometry of an interface can be described by the so-called 'fundamental forms' of the two-dimensional manifold 
/(R). Thus we need to derive the probability distribution of the tangent vectors (i — x,y), the normal vector n 
(see Eq. (j21l), and the metric tensor gij — ti ■ tj. For the tangent vectors 





t, = 1 (47) 



it is sufficient to know the distribution P[dif; a'] of the derivative given by Eq. 1(23 • Notice that we do not normalize 
the tangent vectors because we want to keep the relation gij — ti ■ tj for the metric. Then the probability distribution 
for the metric tensor reads 

1 



P[9xx] = P[9yy\ - , - e (48) 
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FIG. 3: Effective lateral distance a{LKqmax, Lcqmax) at which the variance o"|(a) (Eq. I|48|l ') of height differences equals the 
estimate a'^ag from the local slope variance (see Eqs. (1451 and 13011 ): (a) for aoqmax = 1; (b) for Lcqmax = oo; (c) for ao = 
and (d) for aoqmax ~ 3. The ratio a/ao decreases towards unity for increasing values of Lk and decreasing finite distances ao. 



with g^x = 1 + fxC^) > 1 ^^^d ag — (T3 given by Eq. (|Sn|l . The ofF-diagonal values gxy = gyx = fx{^)fyi'R-) are more 
complicated but can nevertheless be derived explicitly (fx — fo^'^^ and fy = f(^'^\ Eq. H28|l ) 



/oo />oo 
df. / dfyP[fx]P[fy]Sigxy - Ufy) 

(49) 



00 \Jx\ Jx 

00 / 

1 /-dy^-i^^l+i 



27rcr' J y 


1 (\9xy\ 



with (see Eqs.(|551l-(|3IJ) cr| q = Ii,o<^i ~ i""! = ^'^ = '^g/'^- Here, we have used the integral expressions 

dx e-^^-'^/C^-) = J^Kii^) (50) 
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FIG. 4: (a) The probability distribution P[9] of the tilt angle 6 is given by Eq. (I55II . (b) The position 6max of the maximum 
increases with increasing variance = 2a'^ = (Eq. and approaches the limiting value 7r/2 for an infinitely rough 

interface. The maximum value P{9max) depends non-monotonically on ct^. The width A9 at P[0max]/2 reaches its maximum 
value 0.69 for = 0.89. 



and 

oo 

/ - e---'^/^^-) = 2Ko{y%) = -2{K[{y%) + ^^^^) (51) 
J X 



yielding the modified Bessel functions Ko{z) and Ki{z). Then, for the determinant g ~ det{gij) = 1 + (V/(R))^ it 
follows straightforwardly the exponential distribution 



/oo poo 
dh / dfy P[U]P[fy\5{g ~l~fl~ fy) 
-oo J — oo 

1 df. 



"^2.0 "'0 

1 

e "9 



e '"2.0 (52) 



with 5 > 1 and the mean value and variance 

< <? >= 1 + , < .g2 > - < g >2= , (53) 

respectively. Since P[gxy] = ^'[— the mean values of the off-diagonal elements vanish, so that the mean value of 
the metric reads 



_ 1 /l + a'2 
<5^^>-7f(, l + a'2 



(54) 



Using Eqs. (^3), @, and H52|l one can deduce the probability distribution P[9] for the tilt angle 9 between the 
normal n of the interface and the vertical axes e^. With 9 = arccos (^^^ and noting that due to g € [1, oo) one has 
9 e [0, 7r/2), which excludes overhangs, we find (cr^ = 2(t'^ = (t|, Eq. (ISOJ) 

2 sin 61 -i^^ 2 r"'"^ dq ^ksT 
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which is shown in Fig. 2]for various values of a'^. Notice that the maximum of P[6] at 



^max - arctan | ~ 2 + ^da^ ~ Aa^ + 4 



(56) 



approaches the hmiting values 9max = for ag and 9max = 7r/2 for cr^ — s- cx). Close to the vertical direction the 
tilt angle of the normal vector n is distributed linearly in 9 
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1 



P[9^0] = — {l + (-- —)9^ + 0{9 



3 al' 



(57) 



(58) 



whereas one finds an exponentially vanishing distribution towards the horizontal limit 

^ 2^ al{^-9f' 

With the probability distribution given by Eg. J 55 |) one finds straightforwardly the mean value of the tilt angle 9 
(using partial integration and 3.383.10 in Ref. |32|1: 



< 9 > = I dx e ^ arctan \ / xa^ = - — 

2(T„ 



dx X 2 



1 e 



X + an 



V^^(-l)'=2-'=-i(2fc~l)!!af+i , <Jg<l 



k=0 



k=0 



(2fc + l)!! » 



, cr„ > 1 



(59) 



Here, ^{x) denotes the error function 



with its series representations 



and 



TT V^T I TT 2^/n _1 i_ / \ 



X 

2 f ,2 
$(x) = dt 

V'^ J 



2fe^2fe+l 

(2/c + l)!! 



^-(2^ + 1) 



fe=0 



„2fc+l 



fe=0 



(60) 



(61) 



(62) 



with r(fc + 5) = 2-''~(2fc - ^y.l^/^, i.e., r(l/2) = Vtt, r(3/2) = V7r/2, r(5/2) = 30F/4, and T{7/2) = 150F/8. 
Accordingly one finds for the mean squared deviation < (69)'^ >=< 9^ > — < 9 >^ 



< {59Y > = j dx e 1^ arctan \J xafj - < 9 > 



(63) 
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These structural features < 9 > and < (59)'^ > of a fluctuating interface are shown in Fig. [S] as function of the 
variance cr^ = 2a'^. Notice that even in the limit ag oo, i.e., for a soft interface at high temperatures (Eq. H55|l '). 
one does not obtain the mean value and mean squared deviation 



t/2 2tt 



<9>o = I d9 sine I ^ 9 ^ 1 



< {69)'^ = / sin 61 J 




/2 2tt 

2^ 



(64) 



- <9 >l= TT 







of directions n uniformly distributed on the unit semi-sphere. Thus in the limit of a tensionless (i.e., PE{q) = 0) 
fluctuating interface the average defined in Eqs. (|23p - l|25|l is not reduced to the 'random-direction' average < A >q 
due to a geometric coherence induced by the manifold. 




FIG. 5: Mean value < 9 > and mean squared deviation < (56)^ > of the tilt angle between the interface normal n and the 
axis 62 orthogonal to the mean interface (see Fig. Qand Eqs. Q, 15911 . and K'Sl ) as function of the variance of the height 
fluctuations ct^ = ct| = 2a'^ (see Eqs. (gSJ, and The mean squared deviation obtains its maximum value 

< {Sef >ma^^ 0.072 for al = 2.66. 



D. Infrared and ultraviolet limit 



The fluctuation spectrum (see Eqs. l|77|l and ©) 

-^(.) = ^ (65) 

of the interface position /(R) diverges in the limit of large wavelengths, q ^ 0, and vanishing gravitation constant G, 
i.e., Eq ^ Q (see Eq. ^121 ). This infrared divergence is a direct consequence of the occurrence of a so-called Goldstone 
mode with vanishing excitation energy (or effective mass). This mode exists because of the spontaneous symmetry 
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breaking of the translational invariance of the system due to the formation of an interface at /(R) and due to the fact 
that the mean interface position can be shifted without any cost of free energy if G = and without pinning lateral 
boundary conditions. 

In contrast to (y'^{q), the fluctuation spectrum 

-.(9) = .0-^0 ^ (66) 

of the metric (ug — a^] see Eqs. H3Uf) and ^bb\ ). i.e., of the orientations and of the tilt angle, remains finite in the 
limit of a vanishing gravitational constant G (see Eq. (I12|| V It does not diverge because the Goldstone mode at 
q = is flat (corresponding to a rigid shift / = const) and does not contribute to an increase of the fluctuations 
of the local orientations, i. e., fluctuations of the tilt angle 9. In other words, modes with wavevector q contribute 
only proportional to to the distortion of a flat metric which finally suppress the spectrum in Eq. (|66(l in the 
long-wavelength hmit g ^ 0, i.e., a'^{q ^ 0) is a constant. 

In contrast to the infrared case, the ultraviolet limit g — > oo becomes more violent for local orientations as compared 
to the case of positional fluctuations of an interface. Heuristically one may argue that small-scale undulations (ripples) 
do not alter the position considerably but accentuate the more horizontal local orientations due to the many 'steep 
hills' or descents with numerous sharply rising or falling slopes. 

As a consequence the local orientations are very sensitive to the dependence of surface energies at large wavevectors 
- a regime which has been studied recently by grazing incident X-ray scattering [Tol ITlL Il2| . This kind of experiment 
provides only access to o'^(q) which makes it difficult to resolve the details of E{q) at large q. These features are 
enhanced by monitoring the fiuctuations of the tilt angle which are governed by d'^{q) = q^cr'^{q). 



E. Correlation function < 9{0)9{a) > of tilt angles 



The standard characterization of a fluid interface is given in terms of its surface tension, its mean position, the local 
height fiuctuations, and the lateral height-height two-point correlation function. The latter one determines the diffuse 
scattering intensities of X-rays or neutrons from such an interface. In the present study we have put forward the mean 
tilt angle (and its variance) as an alternative characterization of a fiuid (see Fig. |SJ). The temperature dependence of 
the mean tilt angle might be considered to be on a par with that of the surface tension. This alternative description 
in terms of tilt angles reaches the same level of completeness as the aforementioned standard characterization by 
providing the correlation function < 9{0)d{a) > of tilt angles at different lateral positions and a. This goes beyond 
the probability distribution P[0] considered so far. Based on this motivation and inspired by the eventual possibility 
to probe this correlation function also experimentally (see Ref. 22]), in this subsection we determine the correlation 
function < 0{O)9{a.) > and its dependence on the spectrum E{q). It turns out that this correlation function probes 
also non-local properties of the interface. 



|j(n.,™0(a,)} 



has to be calculated, in particular of 
(7a] of the metrics can be 



To this end a family of joint probability distributions P 

the first derivatives f = {fx{0), fy{0), fx{a.), fy{a)). From that the joint distribution P[go 
derived and flnally P[6{0),6{a)] of the tilt angles at a finite distance. 

The concept of orientational fluctuations is frequently used in the context of microemulsions, membranes, and 
macromolcules. There a persistence length is defined which characterizes the correlation < n(0) • n(a) > of two 
normal vectors at finite but not too large distances Is^] . Although related the correlation < 6{0)9{a.) >—< arccos(n(0)- 
e^) arccos(n(a) • e^) > is a more complex function. However, the analytic derivation of the distribution P[9{0), 9{a.)] 
allows us to determine the correlations exactly for arbitrary distances a. 

Following the same procedure as in Eqs. H26() and (|34|l one finds the joint probability of the N values fi := 
/("""'•^(aj) (compare Eq. at the lateral positions a^ (i = 1, . . . , N), 



P 



|jK,m.)(a,)|^^ 



1 



^(27r)^det(S2 



1 



(67) 



with the correlation matrix 



y2 



27r)2 E{q) 

2 



Sni+mi,nj+mj ' (3-i ^j}} — 



(68) 



rii -{-nj /mi+'Tij 



{a.i — SLj) , rii + Tij + rrii + rnj odd , 



17 



S,,^{x) = cos X + ^-^ sin X , 

and the correlation functions for n + m even (compare Eqs. (|29|l and (jSHl) 



(27r)2 E{q) 



(-1) 

oo 

E 

i+i oven 



q"q^ cos(qa) , n + m even , 



(-1)^ 2 



as well as for n + m odd 



'(2^^{qj''''''y 



OO 

E 

i+i odd 



9"^ sin(qa) 
a2(a) 



(-1) 



n+m+i+j + l^x^y 



(69) 



(70) 



(71) 



For the calculation of the correlation function < 9{0)9{a.) > the joint probability distribution is needed only for the 
vector f = (/(i'"'(0),/(''^i)(0),/(i''')(a),/(0'i)(a)) given by the matrices 



y2 _ 



( aio(a) alM\ 

a?:i(a) ^^^(a) 

a|,o(a) aL(a) a'^ 

V^?,i(a) ^0,2 (a) a'2 / 



with the matrix elements 



(72) 



and the determinant 



All 
All 
A22 
Bii 
B12 
B22 



det(S2) 
det(S2) 
det(I]2) 



dct(I]2) 
"det(S2) 



det(I]2) 



(a'4-42(a)-cTti(a)) 
Ki(a)(4o(a)+ 42(a))) 
(a'4-4„(a)-4i(a)) 

(4i(a)a2^2(a) + 4o(a)(a'' - 42(a))) 



(4i(a) (a'4 + 4o(a)42(a)-4i(a))) 
(4i(a)a2_o(a) + al,ie.){a'^ - 4o(a))) 



det(E2) = _ a'4(4_„(a) + 4^,(a)) - 2al,{e.){a'^ + al,{B.)aUe.)) 
+4i(^) +^2,o(a)42(a) • 



(73) 



(74) 
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Repeating the derivation presented in Subsec. IIV Cl one finds the following explicit expression for the joint probability 
distribution 

4singoSinga 1 ff d(f>d(f>' i(taneo.tane.)C(0,0')(taneo,tane,)* fyr. 

"^^'"'"■^J -cos3 0ocos3 0,^d^i/ 47r2 ^^^> 

of the tilt angles at different positions of the random fluctuating interface. The elements of the matrix C(0, (/)') are 
given by 

Cii ((/)) —All cos^ <j) + A22 sin^ + 2Ai2 cos (/) sin 

C22(0') = AiiCos2(/)'+A22sin2 0' + 2Ai2COS(/.'sin0' (76) 
Ci2{(j>, (j)') — C21 = Bii cos cos (/)' + ,622 sin sin cf)' + _Bi2(cos 0sin </>' + sin cos (/>'). 




FIG. 6: Tfie correlation function < 9{0)6{a.) > for the tilt angles 6 depends effectively only on the variance ag of the metric 
and on the functional form of the energy function E{y) = 1 + L^qmaxV^ + '^{y)y^)y i-e., on the wavevector-dependent bending 
rigiditiy K{y). Here, we have chosen (a) = 0, for variances ag = 0.2, 2.0, and 20 as well as (b) ag = 20 and bending 

rigidities K{y) = 1 + {Li^qmax)^y^ with L^qmax ~ 0, 1, and 10 (see Eq. 11841 '). The dependence on the ratio LcQmax of the 
cut-off lengths is small, so that for physical relevant values Lcqmax > 10 no significant change is observed. Here, we have chosen 
Lcqmax ~ 10. The maximum at aqmax ~ does not reflect a minimum in the surface tension 7(g) but seems to be related to the 
sharp cut-off Qmax for capillary waves. 



Finally, the correlation function of the tilt angles 9{0) and 6* (a) reads 

t/2 



< 6l(0)6'(a) > = JJ d0{O)de{a) 9{0)e{a) P[eo,e^] 


(77) 



00 



(ia:i(ia:2XiX2arctan(xi)arctan(x2) _ix*c 



v/d^t(s^ 



-e 



(</,,0')x 



which is shown in Fig. El^a) for the phenomenological capillary wave theory with 7(g) = 70 and in Fig. IHl^b) for 
finite bending rigidity kq (see Eq. (|10|l '). In general, one finds that the centered and normalized correlation function 
(< 9{0)d{a)— < 9 >^)/ < {S9)^ > does not depend sensitively on the particular functional form of the surface 
tension 7(9). In contrast to the one-point probability distribution P[9] given in Eq. IjSSI) . there is only a double- 
integral representation for the two-point probability distribution P[9o, 9a] so that the determination of the correlation 
function < 9{0)9{a) > involves a four-dimensional numerical integration, with the integrand depending parametrically 
on the distance a and functionally on the momentum dependent energy E{q). In order to be able to gain further 
insight into the expression given by Eq. (|77|) we study in the following section the dependence of the relevant variances 
(7^ m(^) (Eq- tZOJ) on the physical parameters entering into E{q). 
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V. DEPENDENCE ON PHYSICAL PARAMETERS 

Within the present Gaussian approach the observables < 6 > and < {66)'^ > are determined by the variance ag 
(see Fig. (S)). In this section we study the dependence of = cr^^g on the physical parameters introduced in Sect. 
mil Since we do not consider systems with finite lateral extensions we set Qmin = and keep G > 0. Expressing the 
momentum dependent surface tension 7(q) in terms of a momentum dependent bending rigidity (see Eqs. ()12f) . H14|l . 
and (Unjl), 

7(g) = 70 + K{q)q^ , (78) 

one has in general 

2 _ I^bT k_i 



(Ji. = 



2^70 



Qmax Sk{Lcqmax] ['«(g)]) (79) 



with the normalized and dimensionless variances 

1 



yk 



Sk{x; [R{q)]) = / dy _ J (80) 







where K{y) = q^^ax^i.^ = <lniaxV)llQ is dimensionless and 



5fc>i(a;^oo;[S = 0]) = . (81) 

A. Phenomenological capillary wave theory 

The phenomenological capillary wave theory is defined by Eqs. (|7jl- (|l()|l with kq = const. For the special case 
Ko = the normalized variances Sk{x, R) (Eqs. (|79|l and H8U|) 'l for fc = 1, 2, and 3 read: 



5i(x;0) = iln(l + a;2) , 

52(x;0) =l_arc^, (g2) 

S,{x;Q) =i(l_l;i(^) . 

For large capillary lengths x = qmaxLc » 1 (Eq. H22|l ). which corresponds to the typical physical situation, one has 
(3.194.5 in Ref. Q) 

TkyiiRo) = Sk>i{x = oo,Ko) = 2^'i(l, -^-TT^; ; -Ko) (83) 

/c — 1 2 2 



where (Eq. 



For k — 2 and fc = 3 one has 



and 



1^0 = ql,ax— = {qmaxL^f . (84) 

70 



72(^0) = 5*2(2; = oo, kq) = ^== arctan v^Ko (85) 



13(^0) = Sz{x = 00, Ko) = ln(l + Ko) , (86) 

2ko 
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respectively. For k = I, Sk diverges in the limit Lc — > oo: 

Si{x; kq) 



1 ^ 2ko(1 + VA) + 1 - A 4ko 
In = , A = 1 — 



2%/A 2ko(1 - %/A) + 1 - A 
^ lnx-iln(l + So) + 0(a;-2) 



(87) 



L« 2 V «09mai: 



with A > 0, i.e., L^, < \Lc. The function 5*1(0;, kq) is shown in Fig. [7| 




La = 

K^max 
kqmax = 1 

T q =2 

kqmax = 5 



5 10 
X = L q 

c" 



20 



c -imax 



FIG. 7: Dependence of the normalized variance crV(fcsT/(27r7o)) = Si{x-Ro = (L„g„„,)2) (Eqs. (|17J, (|HZJ, and l|HHll) on the 
capillary length x — L^qmax in units of the minimum wavelength 2-n /qmax- 



The functions Si{x] kq) and 5*3(2;; kq) determine the variances cr^ = cr^ = = 2cr'^, respectively: 



2 fegJ" c IT - \ 



2vr7o 



and 



with 



2 _ >^bT 2 



27r7o 



^max -J \^ cHmax i 



Ssix; Ko) = ^ Q In [l + (1 + kq)x'^] - Si{x; kq)^ 



(88) 



(89) 



(90) 



The variance (t? characterizes the distribution function of the tilt angle 6 (see Eq. (|55ll and Figs. 0] and |SJ and is 
shown in Fig. |Hl 
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FIG. 8: The dependence of the normahzed variance a'^/{kBTq^g_^/{2-K^o)) — Sz{Lcqmax\ t^o) on the capiUary length Lc = 
^ 70 / Eo (Eqs. I|55|l . H89^ . and l|9()|l ') in units of the minimum wavelength 2Ti/qmax of the capillary waves for various values of 
the length associated with the bending rigidity kq — {L^qmax)^ of the phenomenological theory. S-j, remains finite in the 
limit Lc ^ (X) even for Ro — 0. The variance characterizes the distribution of the tilt angle (Eq. 15511 and Figs. |3]and|SJl. 



B. Microscopic capillary wave theory 

In the previous subsection we have discussed various variances as they are predicted by the phenomenological 
capillary wave theory which is characterized by a constant bending rigidity ko f^nd which leads to a monotonically 
increasing wave length depending surface tension ^{q) = 70 + ko9^ 

The actual microscopic capillary wave theory is defined by Eqs. (|12() . H14I) . and H19() and enters into the variances 
Sk{x; [R,{q)]) via K(q) = (7(9) — 7o9^)/9^ 7^ const. Here 7(9) has a nontrivial structure exhibiting a pronounced 
minimum. In the following we study the influence of this nontrivial structure of 7(9) (or K{q) 7^ const) on the 
variances discussed in Subsec. IV Al and compare this with the predictions of the phenomenological capillary wave 
theory. 

In Fig. 1^ we show the normalized variances 5*1 and S3 corresponding to — and CTj ~ a^, repectively, as 
function of x = Lcqmax (see Eqs. ifTSjl - dSUll V As in the phenomenological case (Figs. [7|and|SJ) diverges ^ Ina; for 
a; — !■ 00 whereas 53 reaches a finite limit T3 — S3{x = 00) for Lc — > 00. For Fig. |51the underlying 7(g) is given by Eq. 
(|19|1 which is characterized by the bulk correlation ^, the molecular diameter tq, and the dimensionless parameter 
Ch', the forms of the functions Sk are independent of 70 = 7(9 = 0) which enters only via Lc (Eq. H22|) 'l. Here and 
in the following we choose Ch = 0.25 which has turned out to describe properly available experimental data of j{q) 
With this choice Sk depends on a; = Lcqmax, roq^ax, ^Qmax, and ^/ro (Eqs. (fT^ and (jSUI)). In Fig. 1^1 we 
choose tq — 2'!T/qmax so that, with £,qmax — 27r<^/ro, Sk{x) depends parametrically only on ^/rp. By using Eq. I|2U|) 
this parametric dependence on S^/vq can be expressed in terms of kq — (27r/ro)^Ko/7o = {L/^qmax)'^ (see Eq. (I84f) ') 
with Kq = "'^f-((C/''o)'' + ^i^/^o)'^) which facilitates the direct comparision with Figs. [7|and|Sl One finds similar 
functional forms but the variances as determined from the microscopic capillary wave theory are larger than their 
phenomenological counterparts. This increase of the variances indicates an increase of interfacial fluctuations due to 
the decrease of the actual surface tension 7(g) for small wavelengths. 

In Fig. IIUI we show the dependence of the limiting value T3 = 6*3 (x = 00) on various microscopic quantities. 
As stated above T3 depends on raq^ax, ^Qmax, and ^/tq where only two of these parameters are independent so 
that Figs. Iior al and llUf b') exhaust the available parameter space. Note, that Fig. [^explores the x dependence 
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FIG. 9: Normalized ~ (a) and S3{x; Lnqmax) (b) corresponding to the variances 

and (jg, respectively, using the microscopic expression for 7(g) (Eq. CHJ) with the choices Ch = 0.25 and ro = 2'K/qmax', 

Ko ~ 1-82 + These variances are larger than their phenomenological counterparts in Figs. [2|and|S| 



within the parameter subspace r^qmax = 27r. Figure ITnT c') shows the dependence on L^^qmax as introduced in Figs. 
□ m and El Using Eq. with Ch = 1/4 one has kq = ^ToJ'o + so that with Eq. |(23) L^q^ax = 

^roqmax ^ + 5 ^ ■ Figure [TUl demonstrates that for Lc 00 the variance for the tilt angle 



as determined from the microscopic capillary wave theory deviates strongly from the predictions of the simple capillary 
wave theory (7(g) = 70) and of the phenomenological capillary wave theory (7(g) = 70 + Koq'^ with kq as given above). 
Moreover, Fig. [TUl shows that depends strongly on the microscopic parameters ^, tq, and qmax- This means that 
the tilt angle distribution is a sensitive probe of the large-q behavior of 7(g) and of molecular details. 

We note that for large correlation lengths the dependence of T3 = S3{x = 00) on ^/ro differs qualitatively for the 
three theories under discussion, because of different functional forms of the bending rigidity function K{y) introduced 
in Eq. (jSUJl: whereas K{y) = and T'3(^) = 1/2 for the simple theory with 7(g) = 70, a constant bending rigidity 

R{y) = Kq leads to T3(^) = ln(l + kq) (see Eq. with kq = ^{roqmax)^ + and therefore to the 

scaling behavior T3(^ >> tq) ^ C *lii(C/''o) for the phenomenological capillary wave theory with 7(g) = 70 + ^og^. 
However, for the full microscopic theory given by Eq. (|19|) one finds for large correlation lengths 

Riy) 0.74Cji^iroqmax)^l-{l + roqrnaxy)e"-»^--y) , (92) 

i.e., K{y) ~ coy^ for <y< raq^^^ ^i*^ ^0 = 0.37C|f (^g^aa;)'^- Thus T3 = = 5^ arct an ^co exhibits 

the scaling behavior T3(^ 00) ^ for the microscopic theory. These differences are important for the temperature 
dependence T ^ which will be discussed in the following section. 



C. Temperature dependence 

The temperature dependence of the local orientations of the liquid-vapor interface is of particular interest. The 
explicit results given above allow one to obtain a transparent view of the analytic structure of < > (T) and 
< ((56*)^ > (T), as well as of its parametric dependence on various features of the intrinsic density profile (such as the 
correlation length ^ and the stiffness parameter Ch) and of the interaction potential w{r) (such as the molecular size 
ro and the range of the potential). This transparency follows from the product approximation for 7(g) (see Eqs. (|14f) 
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FIG. 10: Dependence of the normalized variance T3 — S(x — LcQmax = 00) proportional to the variance ag of the tilt angle 
(Eq. 1911 ') on qmax and on the bulk correlation length ^ (a) and the molecular diameter ro (b) entering the microscopic capillary 
wave theory (Eq. Ch = j. Compa rision is made with the predictions of the simple capillary wave theory (7(g) = 70) 

and the phenomenological capillary wave theory (7(g) = 70 + ^og^ with kq = wloro ( ^ + h^] (Eq. I2()ll 'l: for simplicity 

the phenomenological capillary wave theory is shown only for roqmax = 5 in (a) and for £,/ro = 1 in (b) since curves for other 
values of ro and ^ can be obtained by rescaling the abscissa accordingly. Figure El corresponds to the parameter subspace 

roqmax = 2tv. (c) Plotting Tg as function of L^qmax = \J^roqmax 77 \/l + 5 instead of ^/ro or roqmax allows a comparision 
of the microscopic theory with the phenomenological capillary wave theory for all bending rigidities kq given by a single curve 
(circled solid line). 

and H19|l 'l. We have carried out a full numerical analysis of the temperature dependence for the model considered 
above so that we are able to assess the reliability of the product approximation compared to the full theory described 
inRef. f9|. 

For T ^ Tc the correlation length ^ = CcT^^'^, t = (T^ — T)/Tc, diverges for t ^ where 1/ is a universal bulk 
exponent with v = 0.5 for the present mean field theory (Eqs. IjlSI) and (|16() ). £,q is a nonuniversal amplitude with 
£,0 = ''0/2 for the potential defined by Eq. H18|) . For temperatures well below Tc the correlation lengths and 
^^3-' in the liquid and in the vapor phase, respectively, differ from each other and from the limiting common form 
^ = £,Qt^'^. Although it is straightforward to determine ^'^'^ and ^^^^^ numerically, we have opted for the advantage of 
using the analytic expression given in Eq. ()21|l with — r^jl. The expression in Eq. 121|l has the virtue of fulfilling 
the relations £^^^ > S,{T) > £^^^'> and will be used in the following whenever an explicit expression for the correlation 
length is needed. Moreover, it has the appealing property that at the triple point Tj^ — (2/3)Tc the correlation length 
£,{Ttr) ~ 0.58ro is of the order of the microscopic cutoff length ro, i.e., comparable to the diameter of the particles. 
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FIG. 11: Linear (a) and logarithmic (b) temperature dependences of the variance tTg characterizing the tih angle distribution 
(Eq. I55|l ) for the microscopic capillary wave theory as described in the main text (Eq. lIHSJ); Lcq-max = oo. diverges 
upon approaching the critical point t 0. Within the simple capillary wave theory (7(g) = 70) cTg diverges ~ t~^^'^ whereas 
within the microcscopic theory tTg ~ t"^''^. The crossover to this asymptotic behavior ~ t"^''^ is shifted closer towards Tc for 
increasing values of r^qmax- 



Accordingly, one obtains for the temperature dependence of the variance (Eqs. H79|l . (|20|l . and (|21ll ) 

^9=2A^(l-^j ^3(T) (93) 

with 

3.97 fcsT 2 

A= ^ {qmaxTQ) ■ (94) 

TT'' Wo 

Within the density functional for fluids presented in Sec. Illll one has ^'^'^ — 0.22 so that A ~ Q.Q2d>{qmax'''o)^ ■ 

Whereas for the simple capillary wave theory (7(9) =70) S3{T) = ^ is constant, one finds for the micrsocopic 
theory with the bending rigidity K{y) given by Eq. ll^ the temperature dependence S'3(T — > Tc) ^ t, t = 1 — 

upon approaching the critical point at Tc so that with 70 (t — > 0) ^ t^^^ within mean field theory ag{t — > 0) ^ t^^/^ Int 
diverges for T ^ Te- 
la contrast one finds S^{T Tc) ^ t^\\it for the phenomenological capillary wave theory with 7(9) = 7o + ko9^ and 

kq — ^7o''o + (Eq- (EOl), so that the variance a'^ (Eq. W\\ \ does not diverge but vanishes as ct^ ^ Int 

for T ^ Tc- However, for the microscopic theory one finds S^{T Tc) ^ t which implies ~ i^^/^. 

For the simple capillary wave theory (constant surface energy 7(g) = 70) one has 5*3(7) = i (Eq. (|86|l for = 0) 
so that one finds the temperature dependences (Eqs. (|59|l and H63f) ) 

TT /tT 3 TT 1 6 /tT t / 47r 9 TT 10 [tT 11 , 

< ^ > 2 - V A^' + 2 A^' - V A^' - V 9A^^' + A^" " V A*" + ""^' ^ 

(95) 

<(^^)'>T^ -|^iilnt, t=l-^, 7(g) =70, 

for the mean value < 9 > and the mean squared deviation < {59)"^ > of the tilt angle 9, respectively (see Fig. I12|l . 

The values for the mean tilt angle < 9 > and the mean squared deviations < (69)'^ > depend strongly on the 
microscopic length scales r^qmax, due to the sensitivity on the enhanced fluctuations of the interface on nanometer 
scales as predicted by the microsopic theory. In particular, the actual mean tilt angle is in general larger for the 
microscopic theory than the one predicted by the simple phenomeneological capillary wave theory with a constant 
surface tension 7(5) = 70. Only close to the critical point where the increase of the surface energy 7(9) for large q 
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FIG. 12: Linear ((a), (b)) and logarithmic ((c), (d)) temperature dependences of the mean tilt angle < 9 > [T) (Eq. JSHJ) and 
its mean squared deviation < {SOY > (T) (Eq. 16311 ') for the microscopic capillary wave theory (lines without filled circles) as 
described in the main text and in comparision with the simple capillary wave theory 7(g) = 70 (Eq. I|95|l 'l shown with attached 
circles on the lines for the same values of the cut off length roqmax- The line code is the same for (a)-(d). 



becomes dominant (Eq. (|19|) ). the phenomenological capillary wave theory leads to larger values < 9 > of the mean 
tilt angle, because it neglects bending energies. Due to the decrease in surface energy for wavevectors q ^ 2tt/S^ (see 
Eq. (0) and the increase q for larger wavevectors the asymptotic behavior in the microscopic theory is given by 
7r/2 — < 6' > ~ t"'^/'* and < > ^ t^^^ Ini and not by the expressions given in Eq. (|95|) which hold for the simple 

capillary wave theory. The change of the asymptotic slopes is clearly visible in the double logarithmic plots in Figs. 
IHTc) and (d). 



VI. RELATIONS WITH EXPERIMENTS 



Spectroscopic determinations of molecular orientations at fluid interfaces provide valuable insights into the nature 
and structure of inhomogeneous fluid systems. These techniques render interfacial informations which are comple- 
mentary to those which are obtained by X-ray and neutron scattering 6] which give no direct access to orientational 
ordering at interfaces. Recently Simpson and Rowlen have discussed relationships between the capillary wave spectrum 
and the local tilt angle of the surface 0, Q, 0| . 

One finds in general two different types of approaches to measure local tilt angles of interfaces: 
(i) The interface is decorated with anisotropic particles which keep either a fixed angle with the fluctuating interface 
or the thermal fluctuations of the molecule depend in a systematic way on the interface orientation. Assuming that 
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at low concentrations the decoration does not significantly modify the interfacial fluctuations, one can use the dipole 
moments of these macromolecules as a probe to measure surface specific properties such as the distribution of tilt angles 
(see, for instance, fluorescence and infrared spectroscopy, and Maxwell displacement current (MDC) measurements). 

(ii) The fluctuating interface itself - albeit separating isotropic bulk phases of anisotropic fluid particles (e.g., 
carrying dipole moments) - can exhibit anisotropic features described by orientational profile (see, e.g., Ref. ,341 ) . 
Assuming that such an orientational profile follows rigidly the fluctuations of the interface, one can probe selectively 
the interface, i.e., measure surface specific signals even without adding tracer macromolecules (see, for instance, second 
harmonic generation, AFM micrographs). 



Fluorescence Spectroscopy 



A standard way to measure molecular orientations at interfaces is using fluorescence spectrosopy. For such an 
excitation the polarization of the emitted light depends on the orientation 9,, of the emitting dipole and the orientation 
6a of the absorbing dipole which both are related to the orientation 9 of the molecule at the interface. The relations 
0e{6) and 9a{0) depend on the chemical structure of the fluorescence molecule and are in general fixed and time- 
independent. The so-called time-dependent anisotropy decay is usually defined as r{t) — ^'^~(^«"+^«^)/^ Issl l. where 
the emission intensity lij results from an excitation polarized along the i-axis and measured along the j-axis. Assuming 
that the absorption and emission dipoles are colinear and that the emission intensity is proportional to the projection 
of the dipole onto the detection plane, the measured decay of the anisotropy is given typically by r(t) = (r(0) — 
r(oo))e-*/^ +r(c») (see Refs. [MIMl^ll) and 



— ^4(^y^-2 r(oo) = ,<cos^^.(^)>-,, (96) 



where < • > denotes an average over the orientations 6 of the interface and thus probes the distribution function for 
the tilt angle 9 (Eq. 



Infrared Spectroscopy and linear dichroism 



A classical technique to measure directly the orientation of macromolecules at interfaces is infrared reflection- 
absorption spectroscopy (see, e.g., Ref. (Sql and references therein). Molecular tilt angles of Langmuir-Blodgett 
films can be measured by the comparision of reflection and transmission intensities of infrared beams. The enhance- 
ment factor for the refiection-absorption to transmission-absorption intensities allows for a quantitative evaluation of 
molecular orientations (39[] . However, the dependence on the orientation is rather complex. 

Anisotropic mean orientations of molecules at interfaces give rise to dichroism, i.e., a polarization dependent absorp- 
tion due to the associated orientation of the transition dipole moments. The intensity / of absorption is proportional 
to the mean squared cosine of the angle 9a between the transition moment axis and the electric field: 

/ - < cos^ 9a > ■ (97) 

The angle 9a{9) is determined by the chemical structure of the absorbing molecule and by the orientation 9 of the 
interface normal. Because this relation is fixed, the orientational fluctuations of 9a are given by the probabi lity 
distribution P[9] of the local orientation 9 of the supporting fluid interface (Eq. I55() . For applications see Ref. M 
and references therein. 



Generation of Maxwell displacement current (MDC) 



It is convenient to investigate, for instance, the electrical properties of monolayers of macromolecules floating on 
water surfaces (such as Langmuir-Blodgett films) by the surface potential method with which the Maxwell displacement 
current is measured Ks, 17]. The induced charge Q at the electrode is directly proportional to the vertical component 
niz of the dipole moment of the molecule, i.e., for uniaxial polar molecules 



.7r/2 

Q-<cos6i>= / ^6* sin 61 cos 61 P [61] 
Jo 



(98) 
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Second harmonic generation (SHG) 



As a second-order nonlinear optical process, the frequency doubling of light, i.e., the second harmonic generation, 
is forbidden in isotropic media and thus provides a signal exclusively from spatial regions exhibiting deviations from 
isotropy, such as interfaces. Therefore, it can be used to measure orientational order at interfaces. Such measurements 
have been primarily performed for films of dye molecules of reasonably high symmetry, so that the number of nonzero 
components of the second-order nonlinear optical tensor /J^^^ of the isolated molecule is small. In these cases, with 
one component often dominating the nonlinear response, the nonlinear susceptibility tensor x'^^ of the film can be 
used to extract information on molecular orientation. 

For an uniaxial system of rod-like molecules information about the molecular tilt angle can be inferred from the 
relation (4fl .41. .,42. Al .44 .45,] 

2xizl _ < (sin 0)3 cos > ^gg^ 



(2) ,0^(2) < C0S6' > 



zxx 



The value of each tensor component is proportional to the square root of the second harmonic intensity measured 
under the polarization conditions indicated by the indices. The first subscript denotes the direction of the generated 
polarization driven by polarized electric fields indicated by the second and the third index. 



Direct Measurement of the tilt angle distribution 



Finally, the tilt angle distribution P[0\ (see Eq. (|55|l ) may be measured directly from an AFM micrograph provided 
that the tip diameter is smaller than the typical feature length of the interface. This experiment was done for solid 
interfaces such as silica surfaces (see Refs. JJ5j 14J) but may be achieved also for fluid interfaces when the capillary 
modes can be shock frozen quickly enough. For simple liquids the typical relaxations are too fast for any freezing 
technique to succeed in freezing the capillary waves. But for interfaces involving polymers such as polystyrene- air 
interfaces at temperatures close to the glass transition Tg, the tilt angle distribution may be measured by cooling fast 
below the critical temperature Tg (see Ref. 0|). 



Laser scanning confocal microscopy 



Adding polymers to a colloidal suspension under suitable conditions induces a fluid-fluid dcmixing transition into 
a colloid-rich phase and a colloid-poor phase which are separated by a fluid interface. This phase separation is 
driven by entropic forces. The mesoscopic size of the colloidal particles leads to ultralow values of the surface tension 
which scales up the thermal interface roughness into the regime. The corresponding interface fluctuations can be 
monitored directly in real space by laser scanning confocal microscopy and can be interpreted in terms of the capillary 
wave spectrum |22j| . It seems to be possible to infer from such images also the local interface orientations and thus 
the full tilt angle distribution P{9) as well as the lateral correlation function of tilt angles. This approach would be 
particular promising in that these systems allow one to tune the effective interaction between the colloidal particles 
within a wide range and thus generate designed forms of the wavelength dependent surface tension 7(9) which governs 
the capillary waves. It might be possible to achieve a resolution down to the scale of the colloidal particles, which 
corresponds to short wavelengths on the atomic scale of simple fluids, and thus opens up the possibilities to study the 
influence of the large momentum cutoff. 



VII. SUMMARY AND CONCLUSIONS 



We have studied the local orientations of liquid- vapor interfaces with a particular emphasis on the thermal fluctu- 
ations of the polar tilt angle 9 between the local surface normal and the vertical direction given by the normal of the 
flat mean interface (Fig. Our analysis is based on describing the interface by an effective interface Hamiltonian 
Ti which provides the statistical weight for capillary wavelikc fluctuations. 7i is expressed in terms of a wavelength 
dependent surface tension 7(g = 2tt/\) and we have used three expressions for it: the simple capillary wave theory for 
which 7(5) — 7o is constant and equals the macroscopic surface tension 70 (Eq. (jSJ with kq = 0)j a phenomenological 
expression for ^{q) which takes into account bending rigidities (Eq. H1U|I with kq 7^ 0), and a microscopic expression 
which is derived from density functional theory (Eqs. H14|) and (|19|l 'l. Within the Gaussian approximation we have 
calculated the ensuing variances for the local interface positions and for its derivatives (Eqs. (|26|l and (|30|l 'l as well 



28 



as the variance for the distribution of finite interfacial height differences (Eqs. (|34(l . I|38|l . and (|43|l as well as Figs. 
12 and O))- Based on these considerations we have derived the probability distribution function P[0] for the tilt angle 
(Eq. H55|l ). which exhibits a nonmonotonic dependence on the variance (Fig. ^ and allows us to determine the 
mean tilt angle < 9 > and its mean squared deviation < 9^ > — < 9 >^ (Fig. 13 • In order to probe also non-local 
properties of the interface the correlation function < 9{0)9{a) > of the tilt angles at different lateral positions is 
calculated (see Eq. H77(l ^ which offers an interesting future tool to analyze recent experimental data [22|- The mean 
tilt angle and its moments depend on the physical parameters of the system solely via the variance of the height 
fluctuations ag = <Tk=3 (Eq. (|30|) and Fig. |SJ|. Therefore, Sec. discusses the dependences of the variances (Eq. 
(|79|l 'l on the capillary length Lc (Eq. (|22|l ). the minimum wavelength in/qmax^ the stiffness length (Eq. H22() '). 
and temperature (Figs. I7I11|I . These results demonstrate that these quantities discriminate clearly between the three 
types of wavelength dependent surface tensions mentioned above. This capacity of discrimination between different 
sophistications of effective interface Hamiltonians is transferred also to the temperature dependence of the tilt angle 
< 9 > and its mean squared deviation < 9^ > ~ < 9 >^ (Fig. IT^ . It turns out that the tilt angle is a particularly 
sensitive probe of the wavelength dependent surface tension at short wavelengths which are difficult to reach with 
X-ray scattering. In Sec. IVII we have discussed a variety of experimental techniques which allow one to measure the 
mean tilt angle by optical means. We conclude that the local orientations expressed in terms of the mean tilt angle 
and its moments offer an additional characterization of the thermal states of fluid interfaces on equal footing with 
the surface tension and density profiles and with a particular emphasis on the capillary wavelike fluctuations. This 
underscores the transferability of the concept of the wavelength dependent surface tension from describing interfacial 
structure factors to local orientations of fluid interfaces. 
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